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We evaluate the coordinate space dependence of the matrix elements of the 
commutator of the electromagnetic and gluon currents in the vicinity of the 
light-cone but at large distances within the parton model, DGLAP, the re- 
summation approaches to the small x behaviour of DIS processes, and for the 
Unitarity Bound. We find that an increase of the commutator with relative 
JC ■ distance py as oc (py)f(py,y 2 = t 2 — r 2 ) is the generic property of QCD at 

small but fixed space-time interval y 2 = t 2 — r 2 in perturbative and nonper- 

m 
\o . 

^^ | dipole model (QCD factorization theorem) from the properties of Lorents 

transformation. The increase of f(r) disappeares at central impact parame- 
ters if cross section of DIS may achieve the Unitarity Limit. We argue that 
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such long range forces are hardly consistent with thermodynamic equilibrium 
while a Unitarity Limit may signal equilibration. 



$} , Possible implications of this new long range interaction are briefly dis- 

^ ' cussed. 

> 

X 



Typeset using REVTeX 



*E-mail:blok@physics. technion. ac.il 
TE-mail:fr ankfur@lev.tau.ac.il 



I. INTRODUCTION 

Commutators of the local currents in the coordinate space play an important role in 
the quantum field theory . In particular, they help to visualize the relationship between the 
quantum field theory and statistical mechanics of equilibrium and nonequilibrium systems. 
It has been understood that many properties of deep inelastic processes follow from the 
operator product expansion. In particular, the dependence of the product of local currents: 

< N\j^y)j u (o)\N > = (l/y 2 ) 2 T,PMPy) n < N\O n (0)\N > +NLT terms 

n 

= P^Pu(F(py,y 2 )/(y 2 f) + NLT terms. 

(1.1) 

on the space-time interval y 2 = t 2 — r 2 unambigously follows from this expansion for the 
leading term [1] and from renormalization group. (For certainty we write formulae for the 
product of e.m. or gluon currents and neglect by longitudinal structure functions). We will 
show in this paper that the dependence of operator product on relative distance- (py) at 
fixed y 2 as 

F(py,y 2 ) cc(py)f(py,y 2 ) (1.2) 

with f increasing with (py) follows from the basic properties of QCD. 

The actual behaviour of the structure functions of the nucleon at small Bjorken x = 
Q 2 /2pq is still a challenging question now as it was 30 years ago. (Here — Q 2 is mass 2 
of an incoming photon). So a variety of the new approaches to small x phenomena were 
developed such as: the generalization of the QCD factorization theorem to the amplitudes of 
the hard diffract ive processes which justifies the dipole approach [2,3]; the derivation within 
the BFKL approximation [4] of the dipole approach in the large N c limit, the resummation of 
the pQCD series within the DGLAP approximation [5], the NLO BFKL approximation and 
resummation [6,7], the MacLerran-Venugopalan model [8]; the eikonal approximation where 
a "potential" is evaluated within the DGLAP or BFKL approximations [9-11], Unitarity 
Bound (Black Body Limit) approach [12] 

It is well known that the theoretical description of the high energy processes is sig- 
nificantly simplified in the coordinate space even if actual calculations may appear rather 
cumbersome. The aim of this paper is to evaluate the amplitudes of the deep-inelastic (DIS) 
processes in coordinate space and to visualize the dominant physics. The knowledge of the 
space — time evolution of DIS is especially important for the theoretical description of the 
RHIC program of heavy ion collisions, for the QCD part of the LHC program and for the 
hunt for the new particles [13-15]. 

It has been demonstrated that DGLAP approximation describes well increase of structure 
functions of a proton [16] observed by HI and ZEUS. The experimental data can be fitted 
as: 

xG(x,Q 2 ), F 2 (x,Q 2 ) ocx~ x (1.3) 



with A ~ 0.25 [16]. Basic features of hard diffractive (HT) processes observed by HI and 
ZEUS [16] are well described by the QCD factorization theorem [3]. The success of DGLAP 
at the energy range covered by HERA is due to the energy-momentum conservation law 
restrictions on the possible number of gluons radiated in the multi Regge kinematics. In the 
kinematics covered by HERA this number is equal to 1-2. At even smaller x (this corresponds 
to the LHC kinematics and larger energies) the multiplicity of radiated gluons may achieve 5- 
6 and therefore PQCD approximations become unstable because of the necessity to account 
for large log(xo/x) terms. So various resummation procedures were suggested [5,6]. 

Significant cross section of diffraction in DIS observed by HI and ZEUS [16] are hardly 
consistent with the validity of the leading twist (LT) approximation at large Q 2 but suf- 
ficiently small x. The theoretical analysis shows that the LT approximation is probably 
violated in the kinematics which is not far from that investigated at HERA [16]. Account 
of the conservation of the probability leads to the Unitarity Bound which is the generaliza- 
tion to DIS of the Froissart limit familiar from the hadron-hadron collisions . The Unitarity 
Bound formula shows that the structure functions of a nucleon may increase with the energy 
as [12] 

F 2 , xG oc log 3 (x /a;). (1.4) 

Increase of the structure functions with the energy follows from the increase with the energy 
of the essential impact parameters factor log (xo/x) and from the ultraviolet divergence of 
the renormalization constants in QCD -factor log(xo/x). The conservation of the probability 
permits more moderate increase of the structure functions at the central impact parameters 

as 

F 2 , xG oc \og(x /x). (1.5) 

It has been understood already in 60-s that the dependence of the amplitude of the 
deep-inelastic scattering on v = 2pq and photon virtuality -Q 2 gives unique possibility to 
probe the space-time behaviour of the DIS processes [17]. (On the contrary the amplitudes 
of the soft QCD processes are always on mass shell. So it is impossible to compare with the 
data the dependence of the commutator of currents on the space-time interval y 2 = t 2 — r 2 
and the relative distance py.) Increase with the energy of the coherence length in the target 
rest frame has been suggested in ref. [17], based on the analogy with the QED coherent 
phenomena in the high energy electron interaction in the medium [18,19]. The formulae for 
the coherence length : l c ps Xjlm^x follows from the properties of the Fourier transform in 
[17] since the amplitude of the DIS decreases with increase of Q 2 . Recently, the calculation 
of ref. [21] found that account for the pQCD radiation leads to significantly smaller (but still 
increasing fast with the energy) coherent length as compared to that found within the parton 
model. The first theoretical analysis of the space-time evolution of high energy processes 
in a quantum field theory was given by Gribov [22] and applied for the calculation of the 
nuclear shadowing in the electron -nucleus interactions [23]. 

The dependence of the DIS amplitude on the light-cone interval y 2 was studied exten- 
sively in the parton model [20,24,1,25]. It was also studied in the framework of a Regge 
ansats for small xb structure function behaviour [1,25] . 



The aim of the paper is to show that for for sufficiently small but fixed y 2 the product 
of the currents is the increasing function of py - distance (time), within both the leading log 
(LL) or resummation approaches to pQCD oc (py)f(py,y 2 ) where f is increasing function of 
py at fixed y 2 . Moreover, we shall see that the increase of the product of the color neutral 
currents with (py) is valid for the phenomenological structure functions, describing HERA 
data. Increase of f disappears in the Unitarity Bound approximation but for the fixed impact 
parameters only. In particular, in perturbative QCD 

F(py,y 2 ) oc 9(y 2 )/(y 2 ) 2 (a s N c /^ ggj^ exp(2^a s N c /n\og(py)\og(Q 2 y 2 )) (1.6) 

and in the black limit 

F(py,y 2 ) oc 9(y 2 ) \og 3 ((py))/(y 2 ) 2 ) + peripheral terms (1.7) 

Note that similar increase is characteristic for turbulence. It is well known that the 
velocity-velocity correlator increases with distance in such a system for the case of homo- 
geneous turbulence, i.e. for the scales much smaller than the scale of the entire system. 
[26] . In turbulence such a behaviour arises because the same piece of matter reveals itself in 
different points. Similarly in the deep-inelastic scattering in the target rest frame the same 
dipole reveals itself in different space-time points as a consequence of the large coherence 
length. This explains factor r in the matrix element of commutator. Increase of f{r,y 2 ) 
with r (r ~ t) indicates that produced perturbative system is far from the thermodynamic 
equilibrium. Some caution is to the point: our interest is in the distances less or comparable 
with coherence length. At larger distances deduced formulae are hardly applicable at the 
distances ^> l c where nonperturbative phenomena: confinement of color and phenomenon of 
spontaneously broken chiral symmetry should be important. Discussion of this important 
question is beyond the scope of this publication. Note however that the regime leading to 
the Unitarity Bound corresponds to f not increasing with distance which is a hint for the 
possibility of the equilibrium. 

To visualize physics relevant for the shadowing effects we investigate also the Fourier 
transform of the ratio of the distribution function and the invariant energy s. (Within the 
region of validity of QCD factorization theorem this ratio has a meaning of the cross section 
for the scattering of a dipole off a target [2].) The increase with the distance of a Fourier 
transform of this quantity shows the existence of the long range pQCD interaction between 
the two colorless dipoles. (Remember that the long range interaction related to the zero 
mass of the gluon is cancelled out in the amplitudes of the collisions of color neutral objects 
as the consequence of the gauge invariance). 

Let us note that similar increase with distance can be derived from the calculations of 
structure function evaluated within the Regge pole approximation [1,25] if intercept of Regge 
pole is a(0) > 1 (Only a(0) < 1 case was considered in refs. [1,25]). 

Let us recall conventional definitions concerning the relationship between the products 
of the currents and the structure functions. The structure functions are defined through the 
current product as 

\w, x {q,p) = -(5^ - q li qx/q 2 )W 1 (x, q 2 ) + (l/m 2 )(p M + qj{2x)){p x + q x /{2x))W 2 (1.8) 



W^ x =< p\Jd 4 yex P (iqy)J^y)J x (0)\p > (1.9) 

Here J M is the operator of the electromagnetic current. These structure functions are usually 
redefined into the dimensionless ones: 

F 1 = W 1 , F 2 = {Q 2 /2m 2 x)W 2 (1.10) 

Within the DGLAP approximation these structure functions can be approximated at small 

x as 

F 2 (x,Q 2 ) = f 1 dsG 2 (x/s,Q 2 )g gluon (s,Ql) (1.11) 

J X 

and 

xF x (x,Q 2 ) = I' (ds/s)G 1 (x/s,Q 2 )g gluon (s,Q 2 ) (1.12) 

J X 

The function g , g i UO n(s, <3o) * s the nonperturbative gluon distribution that parametrizes long 
distance contributions while the functions Gi describe the distribution of gluons (sea quarks 
and antiquaries) within the gluon. For the gluon distribution similar convolution formulae 
are valid. (See e.g. ref. [27] for more detailed definitions.) Within the parton model the 
quark -gluon distribution functions are 

G 1 = G 2 = 5(x-s) (1.13) 

For the analysis of the light-cone behaviour it is convenient also to use functions 

V 2 = W 2 /(m 2 N Q 2 ), V L = ((Q 2 /x 2 )W 2 -W 1 )/Q 2 (1.14) 

These functions are free from the kinematic singularities. 

For the theoretical description of high energy processes in the nucleon rest frame it is 
useful to analyze the cross section 

o = F 2 jQ 2 (1.15) 

instead of the IMF parton distribution- F 2 : 

For the gluon-gluon distribution function G this cross-section has the sense of the dipole- 

target cross-section [2]: 

a d = 4ira s x B G/Q 2 (1.16) 

In the framework of the Feynman parton model [20,30,24,1,25] 

V L (x 2 ,px) = -2me(x )5(x 2 )f L (px) 
V 2 (x 2 ,px) = 2irie(x )6(x 2 )f 2 (px) 

(1.17) 



while the calculation based on the Regge models [28] gives 

h ~ (px) a{0) + const, f 2 ~ (px) a(0) - 2 (1.18) 

Our main result is the current-current correlator and the cross-section in the coordinate 
space at fixed and sufficiently small space-time interval y 2 but large relative distances py 
evaluated in QCD using both leading log and resummation models. 

The paper is organized in the following way. In the second chapter we review the results 
of the parton model for the structure functions in coordinate space and show in detail how 
to account for properly the space-time structure of the commutators including causality. In 
the third chapter we evaluate the light- cone correlators of the currents within the DGLAP 
approximation, and find that at fixed space-time interval they increase with the distance 
near the light cone. We found difficult to calculate Fourier transform of amplitude directly 
and to keep causality because of necessity to make approximations. Instead we generalized 
method of calculations developed within parton model in [20,24] . For this aim we found 
convenient to use method of moments including analytic continuation in the vicinity of 
n — > 1. In the fourth chapter we evaluate Fourier transform into coordinate space of the 
phenomenological and theoretical gluon distributions in the small x limit, including both the 
experimental data and the recent resummation models. We also consider the space-time 
behaviour of the structure functions if the Unitarity Bound is achieved at high energies. In 
general we find that the rise of the distributions in the limit x — > leads to the corresponding 
rise of the light-cone product of the local currents. The fifth chapter is the conclusion. 

II. PARTON MODEL IN THE COORDINATE SPACE. 

Let us briefly review the calculations of the structure functions in coordinate space for 
y 2 — ► within the parton model [20,24,1,30]. Within the parton model approximation the 
structure functions are functions of only x. The calculations were carried through in early 
seventies assuming dependence on x as given by Regge formulae, with op(0) <. We need 
to calculate 

f d 4 q 
F(y 2 ,Py) = J -n^F(q 2 ,pq) exp(iqy) (2.1) 

The structure function in the parton model can be derived in the nonperturbative QCD as 
the discontinuity on the cut in the complex pq plane: [24,29]: 

F{q 2 ,pq) = f dx{2pq)e{pq){5{q 2 + 2{pq)x) + 5{q 2 - 2(pq)x)F(x) (2.2) 

Jo 

Here F(x) is the nonperturbative parton distribution in the target. For simplicity we consider 
here spinless quarks. Generalization to spin of quark 1/2 is trivial and does not introduce 
new theoretical phenomena. Let us briefly review the standard way of the calculation of 
these Fourie transforms [24]. We shall start from the integral that is the particular case of 
the integral (2.1)-the integral 



R(y 2 ,py) = / — ^r / dxe(pq)(S(q 2 + 2{pq)x) + <5(g 2 - 2pqx) (2.3) 



(27r)- Jo 

We first calculate the integral 

/(2/ 2 ,M/) = / 7^4e( M )(c% 2 + 2(pg)x) + % 2 - 2^) (2.4) 

This integral is equal to a sum of two integrals, each of them corresponds to the contribution 
of the region pq > 0, and pq < 0: 

I(y 2 , P y) = I + (y 2 ,py) + r(y 2 ,py). (2.5) 

Here 

I ± (x 2 ,px) = / d A q5{q 2 + 2pqx) + 5(g 2 — 2pqx) exp(iqy)9(±qp) (2.6) 

Making substitution q — > q + px in the first term and g — ► q — px in the second we obtain: 

7 ± = 2 cos((M/)a;) / rf 4 g/(27r) 3 5(g 2 )^(±pg) exp(igy) (2.7) 

The latter integral is well known (see e.g. ref. [31]) and immediately gives 

I ± (y 2 ,py) = 2icos(x(py))D ± (y 2 ,py)). (2.8) 

Taking the sum we obtain 

I(y 2 ,PV) = 2icos(x(py))D(y 2 ,py)). (2.9) 

Integrating now over x we obtain the integral (2.3) 

R(y 2 ,py) = 2i dx cos(x (py))D(y 2 ,py) (2.10) 

Jo 

Here the function 

D(y 2 ,py) = (l/(27r))e(H/)(% 2 )) - (mx/(2Vx^))6(x 2 ) J^mxVa?) (2.11) 



is the Pauli- Jordan commutator of the scalar particles (see e.g. ref. [31] for detailed 
analysis of the singular functions in QFT). We retained the full dependence on the nucleon 
mass in order to be sure that there are now singularities in the limit m — ► 0. To account 
spin of quarks one should substitute function D by Green function of spin 1/2 particle -S(y). 
In the LT approximation it is necessary to neglect by masses of quarks. Taking now the 
m — > limit in the latter equation we obtain 

R(y 2 ,Py) = i/ne(py)5(y 2 ) sm(py)/(py)dy. (2.12) 

We can go now to other structure functions discussed in the introduction. The corresponding 
Fourie transforms differs from the integral (2.3) that we had taken by the powers (q 2 ) n (pq) m 



where n and m are integer (but generally non positive numbers). If both n and m are positive 
we can take the relevant integrals just by using the corresponding differential operators. In 
the parton model the scaling leads to the general form of the structure function F(x) and 
one immediately obtains, acting on the eq. (2.11) with the operator pq = —ipd x that near 
the light cone [24] 



F(y 2 ,py) — ► e(py)/(2ir) / dx(Ai(py)5'(y 2 ) cos( x(py)) — 2m 2 xsm(py)5(y^ 

Jo 

+ 2(pd x )(8(y 2 )mx/Jy 2 )Ji(xmJy 2 ))F(x) 



(2.13) 



Here we once again retain nonzero m to be sure there are no singularities. Taking the limit 
m->0we immediately obtain that for the parton model: 



F{y 2 ,py) = (1/tt) / dxe(py) cos(x(py))5'(y 2 )(2(py))F(x) (2.14) 

Jo 

In particular for F(x) = 1 we obtain the Fourie transform 

K(y 2 ,py) = (2/tt) sm((py))5'(y 2 )e(py) (2.15) 

The integral (2.14) is well defined for F(x) ~ x a if a > — 1. This is however not the general 
case. The most interesting structure functions are Fi(x) = l/x,V2(q 2 ,pq) ~ 1/(<7 2 P<?) ~ 
x/q A 1 Vl ~ £><?/<? 4 ~ ^-/( X( l 2 )- It is eas y to see that for these structure function the integral 
(2.14) formally diverges logarithmically and must be regularised. In order to define these 
integrals and satisfy requirement of causality we follow ref . [24] , namely use the differential 
equations: if two functions A and B are connected as 

A = B/q 2 

, then B(q 2 ,pq) = q 2 A(q 2 ,pq), and in coordinate space we obtain 

u A(y 2 ,py) = B(y 2 ,py) 

Let us use this method for the calculation of the structure functions defined above. Let 
us start from iq = 1/x = —2(pq)/q 2 . Then one has 

UF^^px) = +2L(pq) = -2i{pd x )K{x 2 ,px) (2.16) 

Here L means a Fourie image of the corresponding structure function. Since K = 
4sm(py)5'(y 2 )e(py), one obtains the equation 

niq(2/2,M/ = -4i/7r(py)e(py) sm(py)5"(y 2 ) (2.17) 

We look for the solution in the form 

F 1 (y 2 ,py) = A(v)B(u) 



where v = y 2 ,u = py. Since 

DH(v, u) = A{2F V + x 2 F vv + uF uv ) 

where we used p 2 = 0, we obtain 

D(A(v)B(u)) = A{2A v B{u) + vA m B(u) + uB u A v )) 

Note now that 

v S n (v) = -n8 n -\v) (2.18) 

as easily proved by direct calculation. Then if we take A(v) = 5'(v) one obtains 

DH = 4(-5"(v)B(u) + uB u {u)5"{v) 

Comparing this expression with the r.h.s. of eq. (2.17) we obtain 

— i/nu sin(-u) = uB u — B(u) (2-19) 

This differential equation can be easily solved with the boundary condition -B(O) — > 0. The 
solution is 

B(u) = -i/iru / dssm(s)/s (2.20) 

Jo 

Then one obtains 

Fl _+ - l / 7 re(py)5'(y 2 )(py)Si(py) (2.21) 

Here Si is the integral sinus function 

ru 

Si(u) = / sm{s)/sds (2.22) 



(see e.g. ref. [37] for the detailed review of its properties). In the limit of large py that we 
are interested in this paper we obtain 

Fi(y 2 ,py) - -t/2e(py)(py)5'(y 2 ) (2.23) 

Note that the results do not contain any logs as it will follow naively from the corresponding 
diverging integral (2.14) and are causal. 

Exactly in the same way one can calculate V2 and Vl- For V2 one gets 

nV 2 (y 2 ,py) = -i/(27r)e(px)5(x 2 ) sm(py)/(py) (2.24) 

This equation can be easily solved using the ansatz 

V 2 = 6(y 2 )e(py)B(py) (2.25) 

Repeating the steps that lead to the solution of the previous equation we obtain 



uB u + B = sin(«)/« (2.26) 

This equation has the solution 

B(u) = -Si(u) 

u 

We immediately obtain 

V2(y 2 ,py) = -(i)/(27r)(l/py)Si(py)9(y 2 )e(py) (2.27) 



Asymptotically one obtains 



V 2 ~ (-j/4)/(py). (2.28) 



This is just the result of Ioffe [24] (obtained 

practically by the same method). Note that one does not get any large logs using such 
method as one will obtain making naively Fourier transform (see next section). Finally, 
using the same approach one can calculate the function V L ~ (pq)/q 4 - 



The general solution is 



n'V L = 2(py)S (y 2 )sm(py) (2.29) 



V L = (-i/(2n))5(y 2 )e(py) cos(py) (2.30) 



as in ref. [24]. 

Summarizing, in this chapter we reviewed the method of differential equations due to ref. 
[24] of obtaining the Fourie transform of the scaling functions, and stressed that this method 
permits the one to calculate Fourie transforms without violating causality. We considered 
the Fourie transforms of i*\ in the parton model under the condition that i*\ ~ 1/x near 
x — > 0, i.e. has Pomeranchuk behaviour with «p(0) = 1. 

III. COORDINATE SPACE REPRESENTATION OF THE SEA QUARK, GLUON 
DISTRIBUTION FUNCTIONS OF THE GLUON. 

A. Fourie transform of the current-current correlators. 

In the previous section we performed the Fourie transform of the structure functions in 
the parton model. Let us now go to the leading log QCD. Let us make the actual calculation 
for the simplest case- Fourier transform into coordinate space of the structure function of 
the gluons within the gluon. In the case of of quark structure function of a quark or a 
gluon all calculations are practically identical. So there is no need to repeat them. All 
calculations will be made in the target rest frame because the space-time evolution is most 
straightforward in this frame. Our calculations will be legitimate in the limit of the fixed 

space-time interval y 2 but py — > oo. We choose this limit 

because there exists rather direct correspondence between the structure functions at small 
x and Fourier transform. We start from the expression for the gluon structure function which 

9 



is the solution of the DGLAP equation in the double logarithmic approximation [33,34]. To 
derive analytic formulae we neglect the running of the coupling constant which should be 
slow because of the smalness of y 2 The DGLAP equation is 

Q 2 ^G(z, Q 2 ) = (a,/(27r)) j\dz'/z') lGG {z/z')G{z\ Q 2 ). (3.1) 

Here ^qg is the kernel in the QCD evolution equation, 

7 gg = 2N c /z (3.2) 

The solution of this equation is given by 

G(z,Q 2 ) = fdn/(2m)(z ) n ~ 1 /z n (Q 2 /Q 2 r^^ n ' 1 ^ (3.3) 



Here the contour integration over n runs along a straight line parallel to the imaginary axis 
to the right of all singularities of the integral. We use the notation Q 2 = —q 2 if q 2 < and 
Q 2 = q 2 if q 2 > 0. This expression can be rewritten using the scaling variable 

x B = -q 2 /(2pq) (3.4) 

The above solution corresponds to simple initial conditions 

x B G(x B , Q 2 ) = S(x B - x QB ) 

We shall need the Fourie transform of the eq. (3.3). 

G(py, y 2 ) = J d 4 qG(x, Q 2 ) exp(iqy) (3.5) 

Let us first determine the integration area. The structure function is symmetric between 
the u and s channels. So Im s Fi = IhlF u . This means the invariance on the substitution 
x — > —x, pq — > —pq . Thus one can limit the area of integration by the region (pq > 0). In 
this area there are four subregions: 

i. q 2 > and > x^ > — 1 which corresponds to the e + e~ fragmentation into hadrons 
in the field of target. Within the LL approximation this structure function is zero. 

ii. q 2 > and x < — 1 Amplitude in this kinematics can be related with the inclusive 
process :e + e~ — *■ iV + X. We will show that this region gives no significant contribution to 
the kinematics of interest in this paper. In the LL approximation this amplitude is connected 
to DIS amplitude by Gribov— Lipatov relation (see below). 

iii.g 2 < 0, and 1 > x > that corresponds to the deep-inelastic scattering (DIS). 

iv.g 2 < 0, and x > 1. In this area the structure function is because of the energy- 
momentum conservation laws. 

In the second region one has an additional kinematical restriction 

x>l,q 2 >4m 2 x 2 (3.6) 

that just expresses the condition q 2 ^ > q 2 in this kinematical area. 

10 



In the third region it is worthwhile to use instead of q 2 as an invariant variable Q 2 = —q 2 , 
and it is easy to see that kinematically x < 1. 

We shall start from the DIS region. 

Naively, in order to carry the Fourie transform in eq. (3.5), one can use the GIP ap- 
proximation [17,24]. In this approximation one takes into account that the integrand in the 
laboratory reference frame is dominated by 

ql > (g 2 )V(4m 2 ) » \q 2 \ (3.7) 



Correspondingly, one can expand Jq 2 , — q 2 ~ go ~~ <? 2 /(2<7o)- Using this approximation 
one obtains the integrand directly as a function of y 2 . The arising integrals can be easily 
calculated. However, they do not satisfy the causality condition: the commutator is nonzero 
for y 2 < 0, and this condition must be imposed by hand. It is easy to see, taking as an 
pattern the calculations from the previous chapter for the parton model and trying to do 
them explicitly calculating the integrals in the GIP approximation, that the problem is the 
limitation of the integration area by the condition (3.7). Then even when we obtain the 
convergent integrals it is not clear how to obtain casuality naturally. Instead, we shall adopt 
here a different approach. It is possible to prove that the Fourie transform of (q 2 ) n /x m is the 
analytical function of n , uniquely defined by its values in integer n, where the latter function 
is understood as a generalized one. Let us start from the integral (2.15) and multiply the 
integrand by {q 2 ) n - The integral is obtained by acting with the operator D n on the result of 
the integration. For the leading term in the asymptotics in py one obtains: 

tfy _> n»(2/7r)(-ir f dxcos{x{py))8'{y 2 ){py) (3.8) 

Jo 

After differentiating one gets 

q 2n F(x) -> f 1 dx2 2n (py) n+1 S n+1 (y 2 )x n cos(x(py) + mx/2)F{x) (3.9) 

Jo 

There exists however the unique generalized function such that it is an analytic function of 
n and is equal to 5 n (u) for positive integer n [32]. This function is 

J(s,u) = u s - 1 /T(s) (3.10) 

For this function 

lim J(s,u) = 5 n (u) (3.11) 

and we denote u + = 9{u)u (the standard notation in the mathematical literature [32].) Then 
we can extend eq. (3.9) for noninteger n as 

q 2s F(x)/2- 2s -► / dxx- s (py)- s+1 cos(x(px) - sn/2))F(x) 
Jo 

x (4)— 2 /T(-s-l) 

(3.12) 
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Here x\ = 8(x 2 )x 2 . 

Let us shortly discuss the result from the mathematical point of view. It is easy to 
check by using the inverse operator of Laplace as we did in the previous section that the 
asymptotics (3.12) is valid for negative n. Thus we have the problem of restoring function 
that is known for all integer n, analytical in n and has power like asymptotics. Such function 
is uniquely defined [35] , as it is well known from the theory of complex variables (and Regge 
calculus, where the corresponding procedure is called Gribov- Froissart projection). Then eq. 
(3.10) fully defines the function. It is straigtforward to see that this function coincides with 
the one obtained in GIP approximation, except an important difference: we automatically 
achieve causality. Thus our approach-analytically continuing the result of the parton model 
is the only possible approach to the Fourie transform. Once we know how to deal with the 
powers of q 2 , we must put s = a s N c /(ir(n — 1)) and carry the remaining integration over x. 
We obtain using eq. (3.11): 

G{px,x 2 ) = (n/(2py)) f dn/(2m) I dx/x 2 cos(-(a s N c /(n(n - 1)) + 2)vr/2 + xpy)) 
x T(a s Nj(n(n - I)) + 2)(x)^ N ^^ n -^ +2 - n /\(x 2 /(2px))\ asN ^^ n -^ +2 

X (Q 2 \- a sNc/(n(n-l)) z n-l ^2 m -\asN c /(7T(n-l))+2 



(3.13) 



The integration in the above formulae can be easily performed. [37-39]: 

/ sin(a + mxr)x >l ~ 1 dx = (sin(o) (F(imr) + F(—imr)) 
Jo 

— icos(a)(F(imr) — F(—imr)))/(2fi). 

(3.14) 

Here the function F is the confluent hypergeometric function: 

F(x)= l F 1 (fi,/2 + l,x), F(0) = 1 (3.15) 

We are actually interested in the limit of large distances (times), r — > oo. In this limit one 
uses the asymptotics: 

F{imr) = T(fi + 1) exp(i7r/i)/(imr) M + \x exp(irar) / (imr) (3.16) 

F(—imr) = F(fj, + 1) exp(— iirfj,)/(— imr)^ + ^exp(— imr)/(— imr) (3-17) 

or, since 

F(imr) + F(—imr) = 2(T(fi + 1) cos(7r///2)/(mr) M + /isin(mr)/(mr)) (3.18) 

and 

12 



F{imr) — F(—imr) = 2i(T(fi + 1) sin(7r/x/2) /(mr)^ — 2\x cos(mr) / (mr)) (3.19) 

we have 

/ dxx^ 1 sin(a + mxr) ~ 2(r(/i) sin(a + 7r/x/2)/(mr) /i — cos(mr + a)). (3.20) 

Jo 

Using the latter integral we finally obtain the contribution from the deep-inelastic scattering 
region in the Fourie transform. In our case 

a + /xtt/2 = 2a s N c /(w(n - 1))tt/2 - (n + l)vr/2 (3.21) 

H = a s N c /(7r(n - 1)) - n + 1, a = -(+a s N c /(n(n - 1)) + l)vr/2 (3.22) 

Then we have 

G(py,y 2 ) = (ic/(py)) J dn/(2ni)T(a s N c /(n(n - 1)) + 2) 

x (+cos(py - (a s N c /(7v(n - l)))ir/2)/(mr) 

+ T{a s N c /{7i(n - 1)) - n + 1) sin(n7r/2)/ 

/ ((py) asWc/( ^ (n " 1)) "™ +1 |(i/ 2 /(2pi/))r sArc/(7r(n ~ 1))+2 )((5 2 )~ asArc/(7r(n " 1)) 



xzr 1 



(3.23) 

We have found the contribution due to the DIS into the integral over n. The contribution 
due to the fragmentation of e + e~ in the nucleon color field into the integral over n is zero 
in the leading log approximation (see above and ref. [40]). 

Let us now consider the contribution due to the annihilation. The current commutator 
in the annihilation region can be expressed through the correlation function in DIS via the 
Gribov-Lipatov relation [40] 

G a (z,q 2 )=zG(l/z,Q 2 ) (3.24) 

Here z is the Bjorken variable defined in the same way as x, only for the different kinematic 
region. 



G(py,y 2 ) = (ir/(2py))Jdn/(2m)J^ zcos(+(a s N c /(n(n - 1)) + 2) +pyz)) 

x T(a s N c /(7r(n - 1)) + 2 )(z) a ° Nc/Mn - 1)) - n /\(y 2 /2py))\ a ° Nc/ ^ n - 1))+2 

X (Q^-^^/Mn-rD^n-l 

(3.25) 
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The integral over y can be easily taken using the integral 

/ sin(a + mrz)z^~ dz = sin(a + mrz)y^^ 1 dz — / sin(a + mrz)zu~ dz (3.26) 
Ji ' Jo ' Jo 

Then first integral is taken explicitly, while the second was taken above: 

roo rl 

/ sin(a + mrz)z fJ '~ 1 dz = T(/x) sin(/x7r/2 + a) — / sin(a + mrz)z^~ dz (3.27) 

Ji . Jo 

The asymptotic expansion of the latter integral is known, and it is straightforward to see 
that for large r 

oo 

sin(a + mrz)z >l ~ 1 dz ~ + cos(mr + a)/mr (3.28) 

The reason that the latter integral does not depend on mu is that the asymptotics is domi- 
nated by the area z ~ 1, where 

z"' 1 = exp((/i-l)log(^)) ~ 1. 

It is clear that the similar term in eq. (3.13) also comes from the area x ~ 1, and these 
terms correspond to the contribution of the parton model. The two contributions are very 
similar with the only difference that we must use in eq. (3.28) 

a = -(a s N c /(n(n-l)) + 2)n/2 (3.29) 

Then we can write the expression for the structure function that includes both the DIS and 
annihilation regions: 



G(py, y 2 ) = (vr/(mr)) J dn/(2t7r)T(a s N c /(n(n - 1)) + 2)) 

( +{sm(py) cos((a s N c /(n(n - 1)))tt/2) 

+ T(a s N c /(n(n - 1)) + 2)T(a s N c /(n(n - 1)) - n + 1) 

x sin{rni/2)/(py) a ^ /Mn - l)) - n+l 

I \(y 2 /(2py))\ a ° N cl^-^)+^Ql)-^ c /^n-i)) 

(3.30) 

We see that the current-current correlator contains two distinct contributions. The first 
is due to x ~ 1 . The second is solely due to perturbative gluon effects. This part is 
dominated by moderately small x. 

We can now take an integral over n. Let us start from the x ~ 1 contribution 
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(3.31) 



Gi(py,y ) = — (l/(2mr))) / dn/(2ni) cos((a s N c /(ir(n — l)))7r/2) sin(mr)/(mr) 

X r(« s iV c /(7r(n - 1)) + 2)/|t - r |«.tfc/(ir(n-l))+2(g2)-a.JV c /(,r(n-l)) 

x z »- 1 (2m) a ' Arc/(ir(n - 1))+2 

The second contribution is due to the moderately small x It is equal to 
G 2 = T(a s N c /(7i(n - 1)) - n + l)T{a s N c /{Ti{n - 1)) + 2) 

x sin((n)7r/2) 

/ ((mr) asArc//(7r( ™ _1))_n+1 |t - r \ a sN c /(Tr(n-l))+2\^Q2ya s N c /(7T(n-l)) 

(3.32) 

These integrals can be taken using the saddle point method. Consider first the contribution 
where xG ~ const for small x to honor distinctive property of soft QCD amplitudes to sig- 
nificantly more slowly increase with energy as compared to the amplitudes of hard processes. 
This contribution can be studied using the saddle point approximation. Indeed, we have 

Gi{py,y 2 ) = I dn/{2m){l/y 2 ){cos{py) cos{a s N c /{n(n - 1))tt/2) 



( 2 mr/{y 2 QD) cos(« s iV c / (vr(n - l))7r/2)i> s iV c / (vr(n - 1)) + 2)z ) n ~ 1 

(3.33) 
The saddle point is at 

n - 1 = J(a s N/n) \og(2py / y^Q 2 ) / log(z ) (3.34) 

We then obtain: 



G 1 = (cos(py)/(y 2 ) 2 ) cos(yJa 8 N c /7r) log(z )/ \og(2py/y*Q 2 ) 
x T(2 + ^a s N c /n) \og(z )/ \og(2py/y 2 Q 2 )) 
x exp(ya s A^ c /7rlog(mr) hg(z ) \og(y 2 Ql)) 
x (a s N c /^ 4 \og{2py/y 2 Ql) 1 l i / log(^ ) 3/4 



(3.35) 



15 



The last line corresponds to the preexponential. 

Consider now the second integral (3.32). The integral can be rewritten in dimensionless 
variables as 

G 2 = T(a s N c /(7t(n -l))-n+ l)T(a s N c /(n(n - 1)) + 2) sin((n + 1)tt/2)(^)" +1 

/ SQ2 x 2\ as N c /(w(n-l))+2 

x z n - 1 (2) 2 ^ asNc/ ^ n ^ 1))+2) 

(3.36) 
We immediately see that the saddle point is determined from the equation 

\og(pyz ) = -a s N c /(n(n - l) 2 ) log(Q 2 a; 2 ) (3.37) 

We obtain 

n- 1 = ^(a s N c /n)(\og(l/Q 2 x^)/\og(mrz )) (3.38) 

Then substituting in the integral one immediately obtains the asymptotics 



G 2 = r(yJ(a 8 N c /ir(yJlog(pyz )/ log(y 2 Q 2 ) - yJlog(mrzo)/log(y*Q 2 Q )) 

r ( v /(a s iV c /7r)(log(mr)/log(t/ 2 g2) + 2 )) sin( v /« s iV c 7r log(g 2 t/2 )/log( ^ o)/2 )) 

x (a s A^ c /7r) 



a/4 iog(^y 2 ) 1/4 



\og(py) 3 / 4 

x (l/y 2 ) 2 exp(2 v /« 2 iV c /7r log(mr) log(Q§x 2 ))) 

(3.39) 

Let us now check the applicability of the saddle point method. It is easy to see that the 
condition is 

a s log(py)\og(y 2 Q 2 )^l (3.40) 

Thus the equations (3. 37), (3. 38), (3.39) are not valid in the limit a s — > that corresponds 
to parton model. 

We have found the asymptotics of the current-current correlator in the double-logarithmic 
limit of QCD. In this limit the saddle point method is applicable and the correlator increases 
with distances. The applicability condition of this method is evidently the existence of two 
large logarithms: the parameter 

a 8 N c /7rlog(QJh/ 2 )log(py)»l 
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Note, that, generally speaking it is beyond of the accuracy of the method to keep single 
logs in the arguments of the exponents in the above expressilons, and the legitimate answer 
for asymptotics 

G» GLAP (py,y 2 ) = 0(^)^-(a a N e M 1 ^ ] ^^^ eM^a s N c /n\og(py) \og{QW))) 

(3.41) 

Here we put all terms with single log in the arguments of exponents to l.Note that the 
leading asymptotics is given by the integral (3.39). 

Thus we obtained asymptotics for perturbative QCD. Note that delta-function singulari- 
ties on the light cone for F2 and F\ were translated into l/{x 2 + ) 2 behaviour in the perturbative 
QCD. 



B. Coordinate space physics relevant for cross sections . 

In the previous section we discussed the space-time asymptotics of the current-current 
commutator in the LL approximation of PQCD. This commutator has a well defined prob- 
abilistic interpretation in the infinite momentum frame . However in the target rest frame 
significantly more direct interpretation has cross section for the dipole scattering of a target 
and related shadowing effects . The cross-section is equal to the correlator divided by an 
invariant energyi.e. by s, which means the commutator must be multiplied by 2x/Q 2 . Thus, 
in the notations of the introduction we have to calculate 

D(py,y 2 ) = fd 4 qe W (iqy)G(y,Q 2 )/(s) = J d 4 qex P (iqy)(y/Q 2 )G(y,Q 2 ) (3.42) 

Here s = Q 2 /y is the invariant energy squared. 

This quantity is sometimes considered as a "potential" for the interaction between color- 
neutral dipoles. For this quantity we may repeat the analysis of the previous section. It is 
straightforward to see that both for parton model and large Bjorken x, large Q 2 regime of 
DGLAP equations, the effect is the loss of one power of y 2 , in the denominator and the loss 
of one power of (py) in the nominator, i.e. the correlator increases logarithmically in the 
parton model (see previous subsection). For the perturbative QCD asymptotics we obtain: 

D(py,y 2 ) = % 2 )l(a s iV c / 7 r) 1 / 4 l ° i ^ 

(3.43) 
The function D thus increases in the perturbative QCD. 



C. More about parton model. 

The modern definition of the parton model ordinary refers to nonperturbative distribu- 
tions without taking into account perturbative QCD evolution, i.e. for initial conditions for 
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evolution equations, such as DGLAP. The most popular form of the initial conditions (that 
also gives a best agreement with the experimental data) is 

F 2 (x) = C/x a (3.44) 

where a > 0. This case was not considered in early seventies since that time it was assumed 
that a < for the physical cases. The Fourie transform of this function, except the special 
case of the integer a is clearly given by the analytic continuation of the asymptotics obtained 
in the previous section. If we continue analytically the equations from the last section and 
use the results of the section 2 we immediately obtain 

dx cos(x(py))5'(y 2 )(py)x( 1 ~ a ' 1 
o 

~ 2(r(l-a)sin(7ra)G»2/) a 

(3.45) 

Note that for a = this term becomes zero and the asymptotics will be given by 
oscillating term ~ sin(py)/(py) (times the same type of light-cone singularity). 

For the structure function F\ the behaviour is ~ x^ , and one continues analytically 
eq. (3.45) obtaining the increase of the commutator as 

Fi -► 2(r(-a) cos(TTa/2)(py) a+1 5'(y 2 ) (3.46) 

This expression has a pole singularity for a — 0, when we return to the function already 
considered in the framework of the parton model. 

IV. PHENOMENOLOGICAL DISTRIBUTIONS 

A. Small structure functions. 

In the previous section we analyzed the space-time structure of the correlators due to 
a parton model and within the area of applicability of leading order DGLAP equations. 
Including NLO will not change our conclusions. However, at extremly small x (kinematics 
of LHC?) where energy conservation law does not preclude large number of gluon radiations 
in the multiRegge kinematics x limit of the DGLAP equations is, literally speaking, not 
available and instead one needs to use for G either by resummation approaches [5,6], or 
phenomenological ones or phenomenological one from HERA [16]: 

G p (x,Q 2 )~(l/ x r + \Q 2 f (4.1) 

where 

a ~ 0.25, /3~ 0.25 (4.2) 

The theoretical distributions expected within the resummation approaches of [5,6] differ 
from a phenomenological one by logarithmic terms: 



G T (x, Q 2 ) ~ (l/x) 1 - 25 (g 2 ) a25 / V / log(x) 3 (4.3) 

It is straightforward to carry Fourie transform of this distribution. Once again we carry the 
analytical continuation of the parton model formulae. The singularities S(y 2 ) and S'(y 2 ) are 
smoothened into l/d/ 2 .) and l/d/ 2 ^) 2 respectively. We obtain 

G(py,y 2 ) = 7rT(f3 + 2)(pyf +1 /(y 2 ) 2 ^ [' dxx^ " 1 sm(x(py) - fr/2)/(py) 



D(py,y 2 ) = nT((3 + l)(px) /(x 2 ) 1+p f dxx^^ 1 cos(x(py) - [3ir/2)/(py) 

Jo 



(4.4) 

The relevant asymptotics is obtained from the asymptotics of the confluent hypergeometric 
function as in the previous section . For (3 — a < the integrals must be considered in the 
analytic continuation sense, for f3 = a one obtains the logarithmic divergence that must be 
dealt with as in the parton model in the previous subsection. 
Altogether one obtains, if f3 ^ a. 

G(py, y 2 ) = 9(y 2 )7rT(P + 2)F((3 - a) sin(7ra/2))((^) a+1 /{y 2 f +2 (4.5) 

and 

D(py, y 2 ) = 9(y 2 )7rT([3 + 2)T(P - o)(l/2) cos(n(a)/2)((pyr/(y 2 f +1 (4.6) 

If j3 = a 7^ one obtains logarithmic asymptotics: 

G{py, y 2 ) = nF(P + 2) sin(/3vr/2)((^) 1 ^ \og{py) / {y 2 f +2 (4.7) 

and 

J D(^, 2 / 2 )=vrr(/3 + 2)cos(/?7r/2)(^)' 3 log(^)/(t/ 2 )' 3+1 (4.8) 

For the HERA phenomenological case one has (3 ~ a, the same is true for recent phenomeno- 
logical asymptotics due to refs. [5,6]. Thus for them we obtain the logarithmic times power 
increase of D and G functions on the light cone. 

It is interesting to note that recently Ciafaloni et al [6] suggested resummation model 
where the structure function may have a dip in the energy dependence, postponing increase 
to smaller x than in the kinematics of HERA. This will postpone increase of D to larger 
x. Additional factor in the asymptotics ~ 1/log'' (1/x) claimed in ref [6] may change the 



behaviour of D, making it slowly increasing with distance as (py) ' 2 / \/\og(py) and having a 
dip for some interval of py. 

In the case of the leading order BFKL approximation we have (3 = 1/2, a ~ 0.8, and we 
have asymptotics 

G{px,x 2 )~{pyy*/{y 2 f' 2 

D{py,y 2 )~{p y r/{y 2 ?> 2 

(4.9) 
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B. Black limit 

It is worth to analyze coordinate space dependence for Unitarity Limit for structure 
functions for the small x behaviour of structure functions [12] Unitarity Bound= black body 
approximation 

G(x, Q 2 ) oc (l/x)(Q 2 /Ql) log 3 (:r /:r) + peripheral=DGLAP terms (4.10) 

and 

a{s)~ log 3 {x/x Q ) (4.11) 

where x$ is weak function of Q 2 . Doing Fourie transform we obtain 

G(py,y 2 ) oc (py) 9 (y 2 ) log 3 ((py)))/(y 2 ) 3 + peripheral terms (4-12) 

D(py,y 2 ) oc 9 (y 2 ) log 3 ((py))/(y 2 ) 2 ) + peripheral terms (4-13) 

Thus in the black limit both G and D at given impact parameter contain trivial increase 
with distance:factor r is because the same dipole is probed at different space-time points, 
one ln(l/x) is due to ultraviolet divergence of renormalization of e.m. charge. (In (xq/x) 
is due to increase with energy of impact parameters in the scattering process. Excluding 
above factors we find that increase of commutator with distance is stopped within this 
limit. Structure function continues to increase with energy for the configuration in the 
photon wave function with 4k 2 > Q 2 as the consequence of renormalizability of QCD,and 
because of increase with energy of essential impact parameters. 

C. Hard diffractive processes 

Another interesting question is a question of space-time evolution of hard diffractive 
processes. These are the processes 7* + p — > X + p where X is vector meson, dijet etc. For 
this process the relevant amplitude behaves as 

A~(s/Q 2 )l/Q n (4.14) 

where n=l/2 or 1. Repeating calculations as above we obtain that coordinate space ampli- 
tudes increase with distances as amplitudes of LT processes. However dependence on y 2 will 
be weaker by the factor (y 2 ) n . 

V. CONCLUSION. 

We have studied the dependence of the current- current correlators (gluon distributions) 
in coordinate space on py at fixed y 2 close to light-cone. Quite surprisingly, we found 
that all theoretical approaches (DGLAP, BFKL, recent resummation models of small x 
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behaviour [5,6], Unitarity Bound [12] and phenomenological description of data, all lead to 
the increasing with the distance current-current correlators, 

G(py,y 2 )~py/(y 2 ) 2 )(pyr/( y y (5.1) 

here indexes a and (5 are model -dependent at present but positive. The DGLAP equations 
in double logarithmic approximation lead to increase 

~ 9(v 2 )(py)/(v*) 2 (<*.N e /*) 1 ' A ^(^3/4 exp(2 V / a,iV c /7rlog(p ? /) \og(Q 2 y 2 )) (5.2) 

Thus increase the current-current correlators with distance near light-cone is the general fea- 
ture of high-energy scattering processes. Moreover, we see that, apart from the kinematical 
multiplier py this feature appears due to the interaction with the gluons and is absent in the 
parton model, where only fixed (except the kinematical multiplier py )amplitude oscillations 
occur. 

The increase of the commutator with the distance at the light cone is relativistic effect, 
present in the Minkowsky space only, and it is absent in the Euclidean space. 

This feature is closely connected with the known increase of the correlation length at 
high energies (the "Ioffe time", [20,21]). 

Such increase seems to be a characteristic feature of a Pomeron, i.e. of a contribution into 
amplitudes of a t-exchange with vacuum quantum numbers. Fourier transform of amplitude 
with nonvacuum quantum numbers in t channel (contribution into cross section) decreases 
with distance as (l/r) a ^ _1 . Here a is intercept of trajectory of dominant Regge pole 
contribution. 

Let us stress that we consider the asymptotics near the light-cone, i.e. r ~ t. On the 
other hand, for equal time commutator, i.e. small time t but r — > oo dominant contribution 
into Fourier transform arises due to the region of large go ~ 1/i an d small space momenta, 
i.e the region around x=l but q 2 — > oo. In this region Fourier transform oscillates. 

We can also evaluate product of 4 currents in the same way and to obtain similar results 
as above 

j(y)j(z)j(0)j(0) in the kinematics 

y 2 -> 0, z 2 -> 0, r y , r z ->• oo. 

This correlator appears in heavy ion collisions as a correlation function between two hard 
processes which occur at different space-time points. 

Finally, we found that the "potential" function D(y 2 ,py) that has a physical sense of a 
dipole -target "potential" increases with py, in all above cases. This increase disappears for 
the Fourie transform of the Unitarity Bound formula after excluding effects beyond long- 
range dynamics. Thus we have found another instability of the description of the physical 
state in terms of quarks and gluons. The physical consequences of this fact will be discussed 
in more details in the future publications. 
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